Prom laterally modulated two-dimensional electron gas towards artificial graphene 
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Cyclotron resonance has been measured in far-infrared transmission of GaAs/Ala;Gai_a; As het- 
erostructures with an etched hexagonal lateral superlattice. Non-linear dependence of the resonance 
position on magnetic field was observed as well as its splitting into several modes. Our explanation, 
based on a perturbative calculation, describes the observed phenomena as a weak effect of the lateral 
potential on the two-dimensional electron gas. Using this approach, we found a correlation between 
parameters of the lateral patterning and the created effective potential and obtain thus insights on 
how the electronic miniband structure has been tuned. The miniband dispersion was calculated 
using a simplified model and allowed us to formulate four basic criteria that have to be satisfied to 
reach graphene-like physics in such systems. 



PACS numbers: 73.22.Pr, 73.21. Cd, 78.67.Pt 

I. INTRODUCTION 

The range of approaches to explore Dirac fermions in 
condensed-matter physics has recently been extended be- 
yond natural graphenei"— to artificially created lattices 
whose properties such as inter-site coupling or lattice con- 
stant can be tuned. One approach here is to subject gas 
of ultracold atoms to a honeycomb optical lattice^ giving 
rise to Dirac cones in dispersion relations-^^ The same 
has been proposed for lithographically patterned two- 
dimensional electron gases (2DEGs) in semiconductor 
hcterostructures^"— as shown in Fig.[TJ Apart from study- 
ing the Dirac fermions on their own, an appealing per- 
spective in such semiconductor-based systems would be 
to fabricate various proof-of-principle electronic devices 
originally proposed for natural graphene,— for instance 
filters and valveSfii Veselago lenseSfi^ or splittersF^^ The 
electron-beam lithography used to define the artificial 
honeycomb crystal potential, dubbed artificial graphene 
(AG),— allows for much easier control over the device 
details such as the edge geometry, additional or miss- 
ing "atoms" than what would require an atom-by-atom 
manipulatioi>i^ in the case of natural graphene. 

The subject of this paper, fabrication and theoreti- 
cal description of AG, represents a special case of lateral 
semiconductor superlattices (SLs) intensively studied es- 
pecially in the nineties. The body of widely explored 
phenomena can be roughly divided into two classes: 
classical and quantum-mechanical. Magncto-plasmons 
in far-infrared transmission^ and commensurability 
oscillations^^ in a modulated 2DEG can be largely ex- 
plained without invoking quantum mechanics Ji'i^ On 
the other hand, SLs under strong magnetic fields, which 
leads to the magnetic breakdown,— as well as other 
more complex systems2i"— require quantum-mechanical 



ingredients (Bohr-Sommerfeld quantization and tunnel- 
ing between semiclassical orbits). 

In this context, the Dirac cone in the AG spectrum is 
a delicate quantum-mechanical feature just as the Hof- 
stadter butterfl}*2i — Landau bands emerging when p 
quanta of magnetic flux pierce q elementary cells of a 
square SL.— Despite promising recent progress , 
a clear-cut evidence of Dirac fermions in AG has not yet 
been found. Most likely, such an evidence can be pro- 
vided in a magneto-transport or magneto-optical exper- 




FIG. 1. (color online) Artificial graphene: Part (a): hexago- 
nal potential of Eq. ([1]) defining the AG (in units of Vo)- Dark 
regions correspond to carbon atoms in the real graphene lat- 
tice; the primitive cell is indicated, (b): Scanning electron 
microscope image of the surface of one of the samples, (c): its 
atomic layer structure (2DEG depth « 100 nm) and (d): the 
sample layout (not to scale). 
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iment by observing half-integer quantum Hall efFectj^i^i 
or the unique ^Tb scaling in the optical response due to 
inter-Landau level (inter-LL) transitions.—"— More care- 
ful engineering of AG structures is likely needed to realize 
these experiments. 

In this paper, we explore the concept of artificial 
graphene both experimentally and theoretically. We use 
a simplified model of AG to formulate four straightfor- 
ward criteria that must be fulfilled to achieve graphene- 
like physics in conventional 2D semiconductor het- 
erostructures. Among these criteria, we focus on forma- 
tion of a suitable miniband structure, explain how it is 
related to the modulation potential amplitude and show 
how this quantity can be measured in far-infrared mag- 
netotransmission. Although the Dirac fermion physics 
in modulated semiconductor heterostructures has not 
been detected so far (including the experiments presented 
here) , we conclude that simultaneous fulfilment of all the 
four criteria should be technologically feasible. 

II. THEORY 

To arrive at transparent conditions necessary for the 
realization of Dirac fermions in the SL miniband struc- 
ture, we use a maximally simplified AG model which is 
effectively single-parametric The modulation potential 
V{r), r ~ (x, y) shown at Fig. [T{a) is taken as a sum of 
three cosine functions 

y {r) = yo(cos5lr + cosg^r cosg^f), (1) 

where g{ = 2^0(1/73,1), 32 = 2Va(2/%/3, 0), - 
g2 — gl arc the basis vectors in reciprocal space, a is the 
distance between two maxima of V{r) and Vq is the po- 
tential amplitude. Let us note that Vq has to be positive 
to obtain a honeycomb structure rather than a trigonal 
one. The Hamiltonian /2m* + V{r) in basis B of plane 
waves, 

B = ^e'M^.^^y^ = + ^2.92}, (2) 

is a matrix whose diagonal and off-diagonal matrix el- 
ements stand in a ratio determined by Vq, a and the 
electron effective mass m* (in GaAs, 0.067 of the elec- 
tron vacuum mass mg; p is the 2D momentum operator). 
Except for an overall scaling, eigenvalues of the matrix 
depend on a single dimensionless parameter 



This parameter is, up to a factor of the order of unity, 
equal to the ratio between Vo and the kinetic energy Eq 
of a free electron (Vo = 0) at the iiT-point of the Brillouin 
zone. 

Depending on C, we obtain miniband spectra that con- 
tinuously vary from the free 2DEG, through nearly-free 
and more tight-binding-like models, up to nearly fiat 



bands that correspond to practically isolated (artificial) 
atoms. This is illustrated in Fig. [21 where we plot mini- 
band structure for four different values of C,. In Fig. [2^, 
we plot the parabolic dispersion of a free electron (C = 0) 
folded into the newly created Brillouin zone and then fol- 
low the evolution of the miniband structure with increas- 
ing C, namely for C — 0.3, 0.9 and 4. The present Dirac 
cones are marked by vertical arrows and their pseudorel- 
ativistic character has been confirmed by analysing the 
corresponding wavefunctionsi^ 

Importantly, more than one Dirac cone appears within 
the seven lowest lying minibands shown in Fig. [2j This 
fact, to the best of our knowledge so far not mentioned in 
literature, may significantly simplify the quest for pseu- 
dorelativistic physics as we discuss below. The Fermi ve- 
locity corresponding to the lower and upper Dirac cone 
for the modulation strength C, = 0.9 (Fig. [2t) is roughly 
1.4 X lO"' and 3.6 x 10^ m.s~^, respectively. The expected 
Fermi velocity in AG is thus more than order of magni- 
tude lower as compared to natural graphene, where val- 
ues around 10^ m.s~^ are reportedii Another important 
characteristics of the created Dirac cones is their width 
in energy Edc- A closer inspection of Fig. [2] reveals that 
we always get the Eoc smaller than Eq and also that 
Edc strongly depends on C. The typical width of Dirac 
cones, Edc, is dominantly given by two factors: the size 
of the Brillouin zone (i.e. the lattice constant a) and 
the effective mass of the employed semiconductor system 
(which defines the kinetic energy Eq). 

Now we formulate four simple criteria which have to be 
fulfilled to achieve Dirac-like AG physics in hexagonally 
patterned heterostructures. Unless stated otherwise, we 
always consider lattice period a = 200 nm and the ef- 
fective mass of electrons in GaAs m = 0.067mo, which 
match to the samples explored experimentally below. 

(i) Suitable miniband structure: The effective strength 
of the modulation C has to be tuned to get well-separated 
and well-developed Dirac cones. Consistently with our 
calculations in Fig. [51 the range 0.5 < C < 4.0 ensures 
that the Dirac cones do not overlap with other minibands, 
as is the case of low C in Fig.l^b, and also the cones are not 
significantly flattened, which gradually happens for C > 4 
towards the limit of isolated "atoms" . If we strictly limit 
ourselves to the upper Dirac cone, the effective strengths 
( somewhat exceeding 4 are still acceptable. Let us note 
that the chosen range of ( corresponds to the potential 
modulation Vq that ranges from 0.6 to 4.5 nicV. 

(ii) Fermi level positioning and/or carrier density: The 
Fermi level has to be located in the vicinity of a developed 
Dirac cone. If it is well separated from other minibands 
as required by the previous point, the corresponding car- 
rier density is easy to estimate. The number of states 
per unit area in one miniband is equal to 4/(\/3a^) with 
the spin degeneracy included. The carrier density then 
reaches n « 6 x 10^ and 2.5 x 10^° cm~^ with the Fermi 
level located at the Dirac point of the lower and upper 
cone, respectively. If experiments at extremely low den- 
sities (below 10^° cm~^) are to be avoided, we should 
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FIG. 2. (Color online) Minibands generated for several values 
of the parameter (a): = 0, dispersion of a free 2DEG. (b): 

— 0.3, lower Dirac cone develops (indicated by arrow), the 
upper one appears but remains covered by other bands, (c); 

— 0.9, both Dirac cones fully developed, (d): = 4.0, tight- 
binding type narrow minibands form and Dirac cones gradu- 
ally flatten, ultimately becoming again unobservable. Energy 
axis corresponds to GaAs effective mass and a = 200 nm; 
panels a, b, c, and d then correspond to Vb ~ 0.0, 0.4, 1.0 and 
4.5 meV, respectively. 



preferentially focus on the upper Dirac point. Another 
way, if technologically feasible, is to reduce the lattice 
constant. For instance, a = 100 nm implies a very rea- 
sonable carrier density of 1.0 x 10^^ cm~^ for the upper 
cone. As explained above, the miniband structure for 
such reduced lattice constant remains unchanged (after 
rescaling the energy axes in Fig. [3] by a factor of four) , 
provided Vq is increased to keep C constant. 

(iii) Low disorder: The idealized miniband structure, 
as presented in Fig. ^ is in reality smeared out by dis- 
order. The minimal requirement is to have the elec- 
tron mean free path Zg — h^i^J^Trnfe^ significantly ex- 
ceeding the lattice constant a. For 2DEG with den- 
sity of n = 10^^ cm~^ and the relatively low mobil- 
ity 11 = 10^ cmV(V.s), the mean free path reaches 
Zg ~ 500 nm and still remains well above the technologi- 
cally achievable a. We emphasize that this is a necessary 
not a sufficient condition. 

(iv) Cartful probing: A clear evidence for the pres- 
ence of massless particles will likely come from experi- 
ments performed in magnetic fields, using transpor 
or optical^^Ti^ methods. The characteristic spacing of 
LLs, given by the applied magnetic field B and by the 
effective Fermi velocity in AG, then should not exceed 
the Dirac cone width Edc- This condition turns out to 
be numerically close to the requirement that the spac- 
ing of LLs in an unpatterned 2DEG, hwc = heB/m* « 
2 X i3[T] meV, be small in comparison to the modulation 
potential Vg and the width of the particular Dirac cone 
in AG, huJc <C Ejjc- Tolerable magnetic fields are thus 
hundreds of millitesla, since E^c reaches about one meV 



at most in technologically achievable structures. Such a 
low magnetic field requires high quality 2DEG samples. 
If we express this quality in terms of mobility, iJ,B > 1 . It 
is also important to keep temperatures low, kT <^ fkUc, 
which implicates experiments in sub-kelvin range for re- 
alistic AG structures. 



III. EXPERIMENT 

The lateral modulation of 2DEG is in most cases 
achieved either by gating using a specifically patterned 
metallic laye r'^^i^^ or by etching the sample surface using 
methods with high spatial resolutionj ^'^^i^^i^^i'^^ In the 
more common latter case, arrays of dots or antidots are 
fabricatedi^ At a given pitch between dots/antidots, the 
antidot design allows us to get a factor of lower lattice 
constant of AG as compared to array of dots. The etch- 
ing depth serves as a parameter tuning strength of the 
lateral potential. An especially strong modulation can be 
achieved by etching through the 2DEG layer. Such struc- 
tures, with regions fully depleted from electrons?^ have 
been mostly used to study magneto-plasmon effects.— On 
the other hand, quantum effects due to miniband struc- 
ture are typically studied in shallow-etched samples^^ 
such as are subject of our study. 

The studied samples have been prepared by etching 
a shallow array of holes (i.e. antidots) with a triangu- 
lar symmetry. The electron beam lithography and dry 
etching process (Ar+-|- SiC^) have been employed. We 
prepared three samples denoted as A, B and C, see Ta- 
ble m with the etching depth of 15-25, 20 and 48 nm, 
respectively. The diameter of holes was always ss 60 nm 
and hole-to-hole distance, i.e., our AG lattice constant, 
a — 200 nm. Samples A and B have been prepared from a 
wafer with 2DEG in a 20 nm-wide quantum well embed- 
ded between Alo.33Gao.67As barriers and located 100 nm 
below the surface. The electrons in the well are provided 
by two (5-doped Si layers, 15 nm (3 x 10^^ cm^^) and 
25 nm (2 x 10^^ cm~^) deep, see Fig. HJc). The sam- 
ple C was fabricated from a simple GaAs/Alo.ssGao.erAs 
heterojunction located 115 nm below the sample sur- 
face. The triangular well formed at the interface was 
filled by electrons from a Si-doped Alo.33Gao.67As region 
(1.5 x 10"'^* cm""^) separated from the interface by a spacer 
25 nm wide. At all three samples, the lithographically 
patterned area (1x1 mm^) was surrounded by a gold 
frame (50 nm thick) to define the optically active part 
for the transmission experiment, see Fig.[TJd). 

The prepared samples have been studied using the in- 
frared magneto-spectroscopy technique, in this case hav- 
ing a form of the Landau level spectroscopy. To measure 
the magneto-transmission of the sample, the radiation 
of giobar or mercury lamp was modulated by a Fourier 
transform spectrometer. We worked with photon ener- 
gies down to 4 meV at a resolution of 0.125 meV. The 
radiation was delivered via light-pipe optics to the sample 
kept at 2 K inside a superconducting coil and detected 
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by a Si bolometer, placed directly below the sample. All 
measurements were performed in the Faraday configura- 
tion with the magnetic field applied perpendicular to the 
sample layer. 
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FIG. 3. (color online) Relative magneto-transmission spectra 
of studied samples A (dholes ~ 15 — 25 nm, ci2DEG = 100 nm), 
B (dholes = 20 nm, d2DEG = 100 nm) and C (dhoios = 48 nm, 
d2DBG = 115 nm) in panels (a), (b) and (c), respectively. 
The multi-mode character of cyclotron resonance absorption 
vanishes above _B « 5, 6 and 7 T in the sample A, B and C, 
respectively. All spectra are shifted vertically for clarity. 



Characteristic results of the magneto-transmission ex- 
periment are presented in Fig. |21 Relative magneto- 
transmission spectra are shown, i.e. for each photon en- 
ergy, the transmission at a given field normalized to the 
same at i? = T. A well-defined CR absorption of a 
nearly Lorentzian shape is obtained for all three samples 
at higher magnetic fields. When the magnetic field is 
lowered, we observe a more complex behavior with the 
CR absorption split into two or even more modes. The 
energy distance between these modes is clearly different 
in various samples and it is correlated with the depth of 
the etched holes, i.e. with the strength of the modulation 
potential induced by the lateral patterning. An unpat- 
terned reference sample, taken from the wafer used for 
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FIG. 4. (color online) Magneto-transmission spectra (shifted 
vertically for clarity) taken on the sample B at B = 4.5 T with 
three different carrier concentrations adjusted by gradually 
increased illumination time. Splitting of the observed modes 
clearly decreases with increasing carrier density. 



fabrication of the samples A and B, has been also tested. 
As expected, it showed a typical Lorentzian-shaped CR 
absorption at the energy of huc in the whole available 
range of magnetic fields. 

Since the initial carrier concentration at dark was or 
could be modified during technological processing and 
also, since the carrier concentration is affected by the 
near infrared part of the globar/mercury lamp radiation 
(due to persistent photoconductivity), we estimated the 
carrier densities directly from the strength of cyclotron 
resonancei^ This analysis has been performed at higher 
fields, when all specimens provide a well-defined CR ab- 
sorption of Lorentzian shape. The evaluated densities 
for the sample A, B and C are nA,B ^ 0.7 x 10^^ cm~^ 
and nc ^ 1-8 x 10^"'^ cm~^, respectively. Optionally, 
this density could be further increased by illumination 
by an infrared diode. Using the CR absorption width, 
we roughly estimated also the carrier mobility in sam- 
ples after processing which was found to be somewhat in 
excess of 10^ cm^/(V • s) for all three specimens. 



IV. INTERPRETATION 

The observed departure from a single-mode CR at en- 
ergy hcoc shows a clear effect of the lateral modulation 
and we discuss two possible scenarios to interpret this 
finding. The first one relics on a quantum-mechanical 
(single-particle) approach and assumes the AG miniband 
structure due to the lateral periodic potential. These 
minibands are transformed into dispersively broadened 
LLs at non-zero magnetic field. The second scenario, 
a classical one, recalls characteristic multimode CR ab- 
sorption observed in systems dominated by (confined) 
magncto-plasmons. We commence our discussion by 
showing that this latter scenario is, if considered in detail, 
inconsistent with our measurements. 
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Splitting of the cyclotron resonance absorption in di- 
minishing magnetic fields, as shown in Fig. [31 is rem- 
iniscent of magneto-plasma oscillations in a 2DEG.— 
Indeed, the magneto-optical response of our samples, 
which consists from a basic CR absorption line accom- 
panied by one or more modes at (only) higher ener- 
gies, resembles spectra taken on a 2DEG with ID lat- 
eral modulation or on a confined unmodulated 2DEG 
(stripe) Let us therefore compare our results to the 
response of magneto-plasmons in a 2DEG subject to lat- 
eral modulation in both directions, in particular in an- 
tidot latticesJ^i^ In such a case, the far infrared or mi- 
crowave magneto-absorption response includes, among 
others, a characteristic lower branch located below the 
CR energy, which is well documented experimentally 
and explained theoreticallyi^ This lower branch is inter- 
preted as an edge-magneto-plasmon (EMP), circulating 
around the antidot, and its existence is not directly con- 
nected with the symmetry of the antidot lattice.— 

No sign of any EMP mode below CR energy has been 
observed in our experiments. We take this as a clear argu- 
ment against interpretation of the observed multi-mode 
CR response in terms of magneto-plasmons. Another test 
to exclude magneto-plasmon effects was performed using 
external illumination by an infrared diode. Upon grad- 
ual increase of the carrier density n in the sample B (as 
evidenced by the strength of the CR absorption at high 
magnetic fields), the distance of observed modes shown 
in Fig.|4]decreased with n. For magnetoplasmonic excita- 
tions, an opposite trend is expected.— We also recall that 
the experiment of Ref. [l^ concerned deeply etched sam- 
ples and that magnetoplasmon behavior has been typi- 
cally detected at significantly higher densities. 

We now turn to the discussion of experimental facts 
with respect to the four criteria to observe massless 
Dirac fermions in a laterally modulated 2DEG. This pro- 
vides us with basic estimates on how close or far we are 
from Dirac-like conditions. These criteria, formulated in 
Sec. ini are related to the (mini)band electronic struc- 
ture, carrier density, disorder in the system and also to 
the chosen experimental technique and conditions. We 
first discuss criterion (i) for which we need an estimate 
of the modulation potential amplitude Vq. We show that 
infrared magneto-spectroscopy is a method suitable for 
this purpose. 

The main absorption features in Fig. [3] occur close to 
huJc but in contrast to the cyclotron resonance of an un- 
modulated 2DEG, they have an internal structure which 
disappears roughly as 1/B (see Fig. [5]) in the limit of 
high magnetic fields. Such behavior suggests that the 
influence of the modulation potential becomes gradu- 
ally weaker with increasing B and potential energy V{r) 
acts as a perturbation to the LLs. Their spacing, the 
cyclotron energy fiiOc, then provides the dominant en- 
ergy scale compared to the modulation potential Vq and 
the small parameter is Vp/fiwc oc 1/B. In the first- 
order perturbation calculation^ the unperturbed ener- 
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FIG. 5. (color online) Transitions observed in specimens A, 
B and C, cf. Fig. |3] (left vertical axis; depths of etched holes 
are indicated). Vertical arrows indicate magnetic fields at 
which the multi-mode character of CR absorption vanishes. 
The position of the main CR peak in spectra taken on the 
sample A is marked by crosses. The straight line corresponds 
to the theoretical CR-line position with an effective mass of 
m* = 0.067mo, which has been derived from measurements 
on reference (unpatterned) sample (not shown). Lower part 
of the figure (right vertical axis) shows determined values of 
the CR-line splitting AE. The dashed line corresponds to 
the fit of AE for the sample A based on the theoretical model 
discussed in the text. 



gies En = hajc{n + 1/2) become broadened into bands: 

x{2cos/32(^, + ^)cos^+cosH^} 

where = 27r^^g/a^, £q = h/eB and k belongs to the 
hexagonal first magnetic Brillouin zone. Owing to special 
properties of Laguerre polynomials Lnr^ optical tran- 
sition energies, that are En+i,K^,Ky - En,n^,Ky, can be 
rewritten in a simple way. Since we deal with low carrier 
concentrations at which only the lowest LL is occupied, 
we can restrain ourselves to n — 0, 

(5) 

where b{Kx,Ky) denotes the curled bracket of Eq. (U). 

The optical transition energy AEi^q enters the absorp- 
tion probability ai^o for a photon of frequency uj that 
has caused transition between n = and n ~ 1 LLs. 
Provided that the former (latter) LL is full (empty) , this 
probability is proportional to^ 

/ 77r^\{^y'^^^'^y\Px\0,Kx,Hy)\'^S{AEi^o - f^)- (6) 
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Sample 


f-^liolcs 


d2DEG 




c 


A 


15-25 


100 


2.2 meV 


2.4 


B 


20 


100 


3.5 meV 


3.1 


C 


48 


115 


4 meV 


3.6 



TABLE I. Potential amplitude Vb and the corresponding di- 
mensionless parameter ^ for samples A, B and C as derived by 
fitting our data using Eq. (O, see text for details. The etching 
depths and the 2DEG-to-surface distances are also listed. 



If we neglect the dipole transition matrix element (tran- 
sition between the mentioned two LLs is allowed by se- 
lection rules), the characteristic spectral features corre- 
spond to the van Hove singularities indicated in Fig. [51 in 
the occupation- weighted joint density of states (jDOS): 



a(aj) 



5(Ai?i,o-M/o,K..«„(l-/i.«.,.J, 



in which all those transitions at a given energy huj = 
AEifl count where the initial {n = 0) state is occupied 
and the final (n = 1) state is empty, as expressed by the 
Fermi-Dirac factors /. At a filling factor v — nh/eB = 2, 
which was assumed in expression ([6]), a{uj) is a band of 
the width 



w{B) = mfe-'^^'l^ 



(7) 



situated close to — ujc- The width of the band decreases 
with decreasing v (at constant /3), as the filling of the 
n = LL decreases and smaller portions of the magnetic 
Brillouin zone become available for transitions. In the 
limit of very large B, a{yi) turns into a zero-width peak 
at exactly w = 

Eq. ([7]) provides a reasonable basis for interpretation of 
experimental data presented in Figs.[31Ja)-(c). The peak- 
to-peak distance shown as the lower data sets in Fig. [5] 
follows the magnetic field dependence of w(Ii) allowing 
to extract the values of Vb for the particular sample. It 
should be noted however, that the peak splitting observed 
in experiments does not correspond to the full width w 
as calculated using Eq. ([7]) because the lower edge of the 
absorption band is suppressed for v <2. This is the case 
of i? > 2 T and n < 10^^ cm~^ when spin degeneracy 
remains unresolved (Zeeman splitting Sz is roughly an 
order of magnitude smaller than the CR peak width at 
B = 2T;ez= gsehB/2mo « 0.06 meV for \gs\ « 0.44 as 
appropriate in GaAs systems^). Features of the jDOS 
appearing in the absorption band scale as cw{B), where 
< c < 1 is a constant. These features are shown in 
Fig. [S] and correspond to the indicated transitions of the 



broadened Landau bands E„ 



The first states that 



become depopulated upon the filling factor dropping be- 
low two (that is when the magnetic field is increased) 
are those close to the top of the band. Correspondingly, 
the transitions a are the first ones to disappear from the 
absorption spectra. 

For the remaining two features (3 and 7. our form of 
the potential V{x,y) would imply c = 1/9. However, 
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FIG. 6. (Color online) (a) Schematic plot of three lowest ly- 
ing LLs broadened into bands due to the lateral hexagonal 
modulation. Transitions related to the van Hove singularities 
in the joint density of states a(aj) are marked by a, j3 and 
7. Whereas a and 7 are the band edges, /3 is the logarith- 
mic singularity corresponding to the saddle point between two 
neighbouring minima of AiSio in K-space. (b) Magnetic field 
dependence of jDOS a, /3, 7 singularity positions. "Edge" 
shows the lowest in energy allowed transition, as defined by 
the position of the Fermi level (for n — 1.0 x 10^^ cm~^). 
Ec = hwc is the cyclotron energy. 



since the transition /3 gives rise to a logarithmic singu- 
larity which is likely to be smeared out, we focus on an- 
other candidate for an absorption feature: the Fermi edge 
which is also shown in the lower panel of Fig. [5] (transi- 
tions from the states close to Ep to the next Landau 
band). Although the Fermi edge does not precisely scale 
with 'w(B), it appears at frequencies a; « Wc hence we 
take c « 1/3. The values of Vq inferred from fitting our 
data, assuming that the splitting of the CR mode cor- 
responds to ^w{B), are shown in Tab. |T] alongside with 
the corresponding Wc can now return to criterion (i), 
and see that, as anticipated at the beginning of Sec. IIIIl 
shallow etching may create modulation potential favor- 
able for Dirac-fermion physics in AG as quantified by the 
first criterion in Sec. HIl 

While criterion (i) seems satisfied for our samples, the 
carrier density is by a factor « 3 too large compared to 
the requirement (ii), even in the most favourable situ- 
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ation (upper Dirac cone, samples A, B). Since further 
lowering of n may induce metal-to-insulator transitional 
a better strategy seems to reduce the lattice constant. 
Further technological improvements will however, then 
be required to keep criterion (iii) fulfilled. In our case, the 
measured mobility implies mean free path only few times 
larger than the lattice constant and the former is likely 
to deteriorate rapidly upon pushing electron lithography 
closer to its limits of spatial resolution. Once the criteria 
(i-iii) are met for some sample, one should proceed to low 
field {B « 0.1 T) and low temperature (< 1 K) experi- 
ments as required by the last criterion (iv). Since inter- 
LL transitions will then be in sub-THz range and fabri- 
cation of samples homogeneous on a scale comparable to 
the wavelength of absorbed light will be difficult, magne- 
totransport of photoluminescence experiments seem most 
promising. 



SUMMARY AND CONCLUSIONS 

The concept of artificial graphene has been explored 
both experimentally and theoretically. Based on a sim- 
ple theoretical model, we formulated four basic criteria 
that need to be met in order to create and experimentally 
demonstrate the proposed graphene-like bands in modu- 
lated semiconductor hcterostructures. Wc have prepared 
three samples with lateral modulation and investigated 
them using infrared magneto-spectroscopy. The results 
have been discussed with respect to the proposed crite- 
ria with the following conclusions. Etching an antidot 
(hole) lattice on the sample surface creates lateral mod- 
ulation potential with favorable strength, which accord- 
ing to our model, should give rise to a miniband struc- 
ture containing well-developed Dirac cones. An empirical 
rule, which connects the depth of etched holes with the 



strength of the lateral potential, has been found for our 
technological protocol. The quality of prepared speci- 
mens, expressed in terms of mobility or mean free path, 
could be sufficient to resolve the AG electronic struc- 
ture with the Dirac cones, nevertheless, further increase 
of this quality would be desirable. The main obstacle, 
preventing up to now realization of Dirac-like physics in 
2DEG, seems to be related to the interplay between the 
electron density and the lattice constant. Since lowering 
carrier density down to 10^ cm~^ range seems unrealis- 
tic, we instead logically suggest to reduce the lattice con- 
stant below 100 nm. This is technologically challenging, 
given the constraints on mobility, nevertheless still feasi- 
ble. Furthermore, we propose to focus on a higher Dirac 
cone, which is found in AG minibad structure and which 
could be probed at higher carrier densities, namely at 
n « 10^^ cm^^ for a = 100 nm. With these suggestions, 
we believe that massless Dirac fermions can be observed 
in laterally modulated 2DEG, probably using infrared or 
THz magneto-spectroscopy or magneto-transport tech- 
nique, in (possibly not too distant) future. 
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